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In the present work a closed system of kinetic equations is obtained from the truncation of the 
BBGKY hierarchy for the description of the vacuum creation of an electron - positron plasma 
and secondary photons due to a strong laser field. This truncation is performed in the Markovian 
approximation for the one-photon annihilation channel which is accessible due to the presence of the 
strong external field. Estimates of the photon production rate are obtained for different domains of 
laser field parameters (frequency v and field strength E). A huge quantity of optical photons of the 
quasiclassical laser field is necessary to satisfy the conservation laws of the energy and momentum 
of the constituents {e~,e'^ and 7) in this channel. Since the number of these optical photons 
corresponds to the order of perturbation theory, a vanishingly small photon production rate results 
for the optical region and strongly subcritical fields E <^ Ec- In the 7- ray region v < m the 
required number of laser photons is small and the production rate of photons from the one-photon 
annihilation process becomes accessible to observations for subcritical fields E < Ec- In the infrared 
region the photon distribution has a 1/fc spectrum typical for flicker noise. 

PACS numbers: 42.55.Vc, 12.20.-m, 41.60.Cr, 42.55.-f 



I. INTRODUCTION 

The Schwinger effect [l[ of electron-positron pair 
(EPP) production off the vacuum under the action of 
strong electromagnetic fields is one of the few QED ef- 
fects that has not yet been tested experimentally. This is 
due to the huge electric fields E ^ Ec = rn? /e = 1.3-10^^ 
V/cm necessary to observe the EPP production effect in a 
constant external field. Such field strengths are unachiev- 
able for static fields. Therefore, the main attention was 
devoted to the theoretical and experimental study of pair 
creation by time-varying electric fields generated either 
in the focal spot of high-power lasers 0t3 ultrape- 
ripheral heavy- ion collisions [1, 0. Below we will con- 
centrate our attention to EPP creation in strong laser 
fields. The estimations made before in Refs. [H-l^, [lol - [l3j 
showed that pair creation by a single optical laser pulse 
with E Ec could hardly be observed. More optimistic 
results have been obtained for planned X-ray free elec- 
tron lasers (XFEL's) ^3r^ and for counter-propagating 
laser beams in the optical range [T7l - [2l| . 



It is obvious that for subcritical fields E <^ Ec the elec- 
tron - positron excitations have the character of short- 
lived quasiparticles which are not observable after the 
laser signal ceases. In essence, it is a vacuum polariza- 
tion effect. Therefore, the S - matrix methods cannot be 
used [13 and existing estimates [ll] are not reliable. An 
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adequate method is provided by the kinetic theory. Only 
on this basis the different experimental manifestations 
(creation and radiation of annihilation photons 12311 , gen- 
eration of harmonics by focussed laser beams j24| . bire- 
frigency [l^ etc. [2TI. [26H28[ ) find a proper explanation. 
Such a kinetic approach has been introduced in [2^ , for 
a recent review see (soj . However, in these works the 
photon sector was not yet included consistently into the 
quantum kinetic approach. 

In the present work we develop a general kinetic ap- 
proach where both, the Dirac and the photon field are 
consistently included. We base our aproach on the 
BBGKY hierarchy of kinetic equations for the electron 
- positron - photon system generated from the vacuum 
under the action of a time dependent electric field. As 
a first step, we will consider the one-photon annihilation 
process of the quasiparticle EPP, which in the presence 
of a strong field is not forbidden [sij . Estimates of the 
photon production rate for this case constitute the main 
contents of the present work. 

The treatment of the two-photon annihilation channel 
is the next step in the development of the quantum ki- 
netic description on the basis of the BBGKY hierarchy 
for the electron - positron - photon system. The detailed 
treatment of this important class of processes is delegated 
to a subsequent paper. 

Below we will assume an external electric field with 
the 4 - potential (in the Hamiltonian gauge) A'^{t) = 
(0, A(t)) that is spatially homogeneous. It is expected 
that a similar field can be realized experimentally, e.g., 
as a standing wave in the small spatial domain of the 
focal spot of crossed laser beams. The kinetics of quasi- 
particle EPP creation in vacuum has been investigated 
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in detail for the case of a linearly polarized laser field 
within a non-perturbative approach . In the case of a 
rather strong external field Aext(i) some internal field 
Aint(i) will be generated too. The total acting field will 
be equal to A{t) = Aint(i)+Aext(i)j this field is qua- 
siclassical. Fluctuations of the internal electromagnetic 
field on this background can lead to photon excitations 
that, in principle, can be registered outside the active 
zone of the focal spot. 

This article is organized as follows. In Sect. II the sys- 
tem of kinetic equations (KE) in the quasiparticle repre- 
sentation of the electron - positron subsystem in a strong 
time dependent electric field are given. The kinetics of 
the annihilation photons is presented in Sect. Ill and 
the first equation of the BBGKY chain is obtained. In 
Sect. IV the truncation of the BBGKY hierarchy of equa- 
tions on the level of the one-photon annihilation process 
is performed with account of vacuum polarization effects. 
We intend to clearly state all approximation steps and 
therefore indicate them in the text by "Approximation 1" 
... "Approximation 6" in brackets at the place of their 
introduction. The spectrum of annihilation photons is 
investigated here for the large (multiphoton processes) 
and small (strong fields) values of the adiabaticity pa- 
rameter. Estimates for the photon radiation rate show 
that it is very small in the case of optical lasers with 
subcritical fields E <^ Ec since a huge number of optical 
photons from the quasiclassical laser field is necessary 
here to overcome the energy gap. However, the effect 
becomes quite observable in the domain of strong fields 
E 5 . Ec and high frequencies (e.g., for the XFEL domain 
[14| ) . Finally, we summarize the basic results of this work 
in Sect. V. 

We use the metric g^i, ~ diag(I, —I, —I, —1) and nat- 
ural units h = c = 1. 



II. KINETIC EQUATIONS IN THE ELECTRON 
- POSITRON SECTOR 



In general, the complete system of equations for a con- 
sistent kinetic description of the electron - positron - pho- 
ton plasma consists of 

• the KE's for the distribution functions fe,p{p,t) of 
the electron and positron quasiparticle components 
in the presence of the strong quasiclassical total 
Maxwell field A{t) = {0,0, A{t)); 

• the KE for the distribution function of the photon 
component, and 

• the Maxwell equation for the quasiclassical internal 
field Ai„t(i) . 



We start from the standard QED Lagrangian C = Cqc + 
Cj taking into account the interaction of the electron- 
positron Dirac spinor fields with a quasiclassical electro- 
magnetic field in Cqc and with a quantized one in Cj, 



{iYD^^ - (Z?»7^V'} - m^V', (2) 



(3) 



where D,, 



ieA„ 



'n — I -^-^ •-fj.it) ■ Thus, the spatially homo- 
geneous nonstationary quasiclassical Maxwell field A{t) 
is a background to the quantized photon field A^{x) = 
(0, A(a;)). We assume that the intensity of the quantized 
field is rather weak and its influence on the state of the 
system can be neglected. In other words, the electron- 
positron system plays the role of a photon source. The 
photon in - vacuum \in) is defined such that for this state 

{A^ix))=o.^ 

Below we will not consider the backrcaction problem 
because for subcritical fields E ^ Ec the internal field is 
negligible, and hence E{t) = Ecxt{t)- 

Thus, the idea of the work is the following: With the 
fermion KE's (see below), in principle, the quasiparticle 
distribution function ([T]) can be determined, taking into 
account the spin degrees of freedom. The separation of 
the quasiclassical Maxwell field allows to formulate a non- 
perturbative approach to the photon kinetics (Sect. III). 
In the present work, the photon spectral distribution will 
be considered mainly by including the one-photon anni- 
hilation process only (Sect. IV). 

The kinetics of electron - positron vacuum pair creation 
under the action of a linearly polarized electric field has 
been studied in several works, (see, e.g., [l, 32, 33] and 
the references therein). The corresponding generalization 
to the case of an arbitrarily polarized time-dependent 
electric field was obtained in j34l - l36j . In the present work, 
we will use the oscillator representation ^ [32| which leads 
to the nonstationary spinor basis [33 - [36j 



B(p)[lu+,0,P^P.] , 
i?(p)[0,c^+,P+,-p3] , 

B{p)[^P\-P^,u;+,0] , 
t{-p,t) = B{p)[-P+,p-\0,u;+] 



u+{p,t) 
vt{-p,t) 



(4) 



with the usual orthonormalizion conditions j42| 
ut{p,t)vi3i~p,t) = , 

ut{P,thp{p,t) = V+{-p,t)vp{~p,t) = Sap , 



Ua{p,t)ui3{p,t) 
Va{p,t)vi3{p,t) = 



uj{p,t) 
m 

u}{p,t) 



Sa 



/3 



(5) 



We assume that the electro-neutrality condition holds 

fe{p,t) ^ fpi^p,t) ^ f{p,t) . (1) 



^ In this representation the diagonalization of the Hamiltonian in 
the presence of an external field is achieved at once, without the 
canonical Bogoliubov transformation procedure. 
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where Ld{p,t) = y/r. 



p - eA(t), P± 



= Lo + m and B{p) = [2a;w+]^^/^. The 
Hamiltonian has the diagonal form in this quasiparticle 
representation 

-b^{-p,t)b+{-p,t)] . (6) 

Finally, the Dirac equation in the presence of the exter- 
nal quasiparticle field Acxt(0 generates the Heisenberg - 
like equations of motion for the time-dependent creation 
and annihilation operators 

a{p,t) = -[/(i)(p,0a(p,i)-f/(2)(P,i)&+(-P,i) 
-iuj{p,t)a{p,t), 

b{-P,t) = 6(-p,t)C/(i)(p,t) + a+(p,0t/(2)(P,i) 

-zw(p,i)6(-p,i) , (7) 

with the matrices in the representation @ being 

[/(l)(p,t) = iwC(p)[PE]cr = iUcr, (8) 
C/(2)(p,<) = q<T = C(p)[P(PE) - Elulu+](t . (9) 

Here, E(t) = —A{t) is the strength of the external elec- 
tric field for |Aint| < |Aext| and C(p) = e/{2uj'^uj+). The 
matrices (O and © describe the different vacuum effects 
in the presence of an external electric field (polarization, 
spin rotation, EPP creation). 

The KE's for the electron - positron component of the 
plasma follow from the equations of motion ((T]) and the 
definitions of the electron and positron distribution func- 
tions in the instantaneous representation 



faf3iP,t) = (a+(p,t)a„(p,t)) , 

fapiP,t) = {bpi-p,t)bti-p,t)) , 

and also two additional functions 

/l?(p,t) = (4(p,t)6i(-p,t)), 

fi'pHp^t) = {b0{-p,t)aa{p,t)} , 



(10) 



(11) 



describing vacuum polarization. The KE system in ma- 
trix notation is then [SJ-ISal 



r = [r,c/(i)] + (/(+)c/(2) + c/(2)/<-') , 

/^+^ = [/^+^C/(l)] + (C/(2)/-rC/(2))+2^^./^+^ 

f^-^ = [f^-\U^i)] + (/C/(2) - t/(2).r) - 2zc./(-) . 

(12) 

If the standard decomposition in the basis of Pauli ma- 
trices is used (fc = 1,2,3) 



where /o = ^Tr/ , A- = ^Tr/crfc, the KE's (HH) can be 
rewritten in the spin representation, where the first set 
of equations concerns the spinless distributions 

/o = -2qu, 

/o = 2qu , 

uo = 2wuo + (f - f'')q , 

vq = -2ujuo , (14) 

and the KE's for the spin distribution functions arc col- 
lected in the second set 



fk 

tk 

Uk 
Vk 



-2uo%-2[fU]fc + 2[vq]fe , 
2wogfe-2[rU]fe-2[vq]fc , 
2L0Vk - 2[uU]fc + (/o - f^)qk , 
-2c^^ife-2[vU]fe + [(f + r)q]fc 



(15) 



Some simple applications of this system of KE's can be 
found in [34|-|36|. 

In the simplest case of the linearly polarized photon 
field A.{t) — (0,0, A(i)) the system of equations is 
transformed into a coupled system of three ordinary first- 
order differential equations [2^ [30| 

1 



/ = ' 
w A(l - 2/) - 2lov 
V = 2lou , 



(16) 



where X{p,t) = eE{t)e±/uj'^, £±ip) = + p"^ and 

p^ = (p'^\p(^', 0). The system (|16p corresponds to a 
KE of the non-Markovian type 

/(p, t) = iA(p, t) dt'A(p, t')[l - 2/(p, ^')] 



2^ d™(p,T) 



(17) 



This KE and its representation (\16^ have been used in 
numerous applications, see [30j . 

The further development of the kinetic theory ap- 
proach to vacuum particle creation is progressing in 
search of different exact and approximate solutions of 
the problem (see, e.g., fSTl - lioj) and also in extending the 



possibilities of the formalism 



) an 



/ = /O + /fcCTfe 



(13) 



III. PHOTON SECTOR 

A. The basic equations 

An external electric field that generates an instability 
of the vacuum with respect to electron - positron pair cre- 
ation is accompanied by the appearance of internal cur- 
rents and electromagnetic fields (back reaction problem) . 
Quantum fluctuations of this internal flcld are interpreted 
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as photon excitations, which can leave the active zone 
(focal spot) and can be registered experimentally. Below 
the KE for the description of the photon component will 
be obtained and investigated for some simple situations. 

For the construction of the photon kinetics it is neces- 
sary to develop the corresponding generalization of the 



J 



quasiparticlc formalism developed in Sect. II using it as 
a nonperturbative basis. 

Let us introduce at first the interaction with the quan- 
tized electromagnetic field (x) in the fermion sector of 
the theory by means of the substitution Hf — > Hj + Hiat 
in the Heisenberg-like equation of motion ([7]) 



a(p,t) + C/(i)(p,t)a(p,t) + t/(2)(p,t)&+(-p,t) 



-i[a{p,t),Hf + Hi 



6(-p, t) - b{-p, t)[/(i)(p, t) + a+(p, t)C/(2) (p, t) = -i[b{-p, t), Hf + i/i„t] 

I 



(18) 



where Hf is the Hamiltonian ^ of the fermion field 
in the quasiparticle representation and Hmt is the usual 
Hamiltonian of the interaction with the quantized field 



i^int(i) d^x: ^ix)-f''A^{x)^{x) 



(19) 



The time dependence of Hi^t {t) is due to the nonstation- 
arity of the system that is reflected in the characteristic 
oscillator representation of the decomposition of the field 
operators ■(/'(x), ip~^{x) in the nonstationary basis ([4]), ([5]). 
The same source (external field) generates the nonsta- 
tionarity of the quantized electromagnetic field. How- 
ever, this does not induce an alteration of the mass shell 
of the photon field, fc^/c^ (in contrast to the electron 
- positron field, w(k,t)), and then the standard decom- 
position is valid. 



A^{x) = (2^)-3/2 



— zkx 



(20) 



where A^(k,t) = A'-^\k,t) + Al ^(-k,i) with the condi- 
tion {^^,^^(k,t)}+ = A''^\-'kJ) and the standard com- 
mutation relations 



[4-)(k,t),4+)(k',t)] = -5^,<5(k-k') 



(21) 

The Hamiltonian of the interaction (jl9p in the oscilla- 
tor representation has the form 

ffint(0 = e(2^)-3/2^ J d'p,d'p,^ 

^{Pi - P2 + k) 
: { M^a (Pi , P2 , k; t)a+ (pi , t)ap (p2 , t) 
+ M/ja(Pi> P2, k; t)a^{pi,t)b+{-p2,t) 

+ M Pa (Pl , P2 , k; t)6a (-pi , t)af3 (P2 , t) 

+ M^a(Pi> P2, k; t)bai-pi,t)b+{-p2,t)} 
Ar.{k,t) : . (22) 



correspond to the momenta of photons. The spinor con- 
volutions are 



N]^a(Pi,P2,k;t)=^„(pi,t)7%(p2,t)e:(k) , (23) 



where the spinors (p, t) and rja (p, t) are taken from the 
set {mq(p, i), Uq(— p, t)} (jU, and e^,e^ are the photon 
polarization unit vectors, while e = k/A: and — 6'^. 
The relationship [HI 



4±)(k,i)-E^M4^^('''*) 



(24) 



r=0 



has been used here. 

The Hamiltonian of the free photon field is given by 



Hph {t) 



^ f d^k kAi+^ (k, t)Al^^ (k, t). 

^ — 1 n ^ 



(25) 



Here and below the vectors pi, p2, . . . are used for nota- 
tion of the canonical momenta of fermions and ki , k2 , . . . 



The system of Heisenberg-like equations of motion 
with account of the photon subsystem (jlSp can be writ- 
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ten now in explicit form. For example, This is the exact system of equations with nonpcrturba- 

tive account of the external electric field. 

^[/(2)(p,t)5+(-p,t) 

f (Pk 

+ie(27r)^''^ / d^pi—= S{p pi + k) B. The first equation of tiie BBGKY iiierarchy 

J V 2fc 

|a(pi,i)[uu]'"(p,pi,k;i) -y^^ ^^^^^.^ ^^^^^ ^^^e single - time, two - point photon 

,\r- ir/ 1 ,\1 ^ /, ,\ correlation function in momentmn space 

+b{-pi,t)[uvY{p,Puk;t)jAr{k,t) , (26) 

b{~p,t) = -iuj{p,t)b{-p,t) + b{-p,t)U(i-,{p,t) Frr'iKk'.t) ^ {Al+Hk,t)Al7\k',t)) . (29) 

h2)i 



+a+(p,t)L/(2)(p,t) 



+ [vv 



3/2 f 3 ^ spatially homogeneous case it is diagonal m the 
-*e(2^)-^/^ 

%/2fc ^^^^ P + momentum k and polarization r with the photon distri- 

(- bution function Fr{k,t) as matrix element, 
\^[uvY {pi,p,k;t)a+ {pi,t) 

npi,p,k;t)b{~pi,t)}Ar{k,t) , (27) {Ai+\k,t)Al7\k',t))=Srr'S{k-k')Fr{k,t) . (30) 



The diagonalization in the polarization indices is an ap- 

Te(27r)-3/2^_ f ^-p^d^p^ j(p^ _ p2 ip k) proximation here (Approximation 1). 

"^/^fc J Then from the commutation relation (PT|) it follows 
|a+ (pi, t)[mt]''(pi, P2, k; t)a{p2,t) 

+a+(pi,t)H'^(pi,P2,k;t)6+(-p2,i) (4-Hk,04^^(k',t)) = 5...<5(k-k'){l+F.(k, t)}. (31) 
+fe(-pi, t)[wu]''(pi, P2, k; t)a(p2, t) 

,\r- ir/ 1 + / Let us write now the first equation of the BBGKY hi- 

+o(— Pi, (pi, P2, k; t)o^(— P2, t) > . , r i , . „ . jt^ttty. 

J erarchy tor the correlation function (|29p using the photon 

(28) equation of motion ([28]) 



J 



Frr'ik,k',t) = ^e(2^)-3/2V fd^Pld^P2\ - ^(5(pi - P2 - k) 



3a (Pi , P2, k; (a+ (pi , t)a^ (p2 , t)Al, ^ (k', i)) 



+ [uv]l^{pi,P2,k;t){a+{put)b+{-p2,t)Al, \k',t)) + [H^«(Pi, P2, k; t)(6„(-pi, t)a^(p2, '(k',0) 
+ H0a(Pi, P2, k; t){b^{~pi,t)b+{~p2,t)Al7'> (k', t)) 



^ =5(pi-p2+k') 



'2k' 

[uuY^^{pi,P2X ■,t){at {pi, t)af,{p2,t)Ai+Hk,t)) + [uvYp^^^ 
+ M^'a(Pi> P2, k'; t)(&„(-pi, i)«MP2, t)4+^ (k, t)) + H^'jPi, P2, k'; t){b,,i-p,,t)b+{-p2,t)Ai+\k, t)) | 

+i{k- k')Frr'{k,k',t) . (32) 



The last term on the r.h.s. of Eq. (j32[) describes the quan- forbidden in absence of an external field become possible 

tum beating of two- photon states and can be omitted in here: e.g., in the lowest order of the perturbation thc- 

the approximation (j30p . ory it is the one-photon annihilation, the simultaneous 

Thus, the kinetics of the photon states is defined by creation of an electron - positron pair and a photon 

the different forced processes of either one-photon scat- Apparently, the latter processes are strongly suppressed 
tering of electron and positron (considered as quasiparti- ^^e region of subcritical fields. In highest orders of 

cles) or their creation and annihilation. Some processes Perturbation theory the number of this kind of forbid- 
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den processes increases abruptly. In the general case, the 
necessity of taking into account such kind of processes 
strongly complicates the problem. 

According to the approach given in Sect. II, the in- 
teraction of the photon subsystem with the electron- 
positron one can be taken into account in the frame- 
work of the standard perturbation theory with the fine- 
structure constant a = e^/ (47r) as an expansion parame- 
ter. This justifies the truncation of the BBGKY hierar- 
chy for the correlators occurring in Eq. ()32p . 



For obtaining a closed photon KE it is necessary to 
perform some truncation procedure for the correlators 
occurring in this equation. The simplest truncation of 
the type 



IV. ONE-PHOTON ANNIHILATION CHANNEL 

A. Annihilation channel 

Below we shall restrict the discussion to the annihila- 
tion channel only. Keeping the relevant terms in Eq. ((32|) 
('incoming - outgoing' terms, the second line in the first 
square bracket in r.h.s. of Eq. (|32p and the third line in 
the second bracket), we obtain 



(6„(-pi,t)a/3(p2,t)4^^(k,t)) ~ 

(6„(-pi,t)a^(p2,t))(4±)(k',i)> =0 (34) 



Fj-j-' (k, k , t) 



{^='5(Pi - P2 - k)[uvYp^{pi,p2,k;t) 



{b+i-p2,t)a+iput)47\k',t)) 

-i=(5(pi - p2 + k')[vuY^^{pi,p2,k';t) 
V 2fc 

(fe„(-pi,t)a0(p2,t)4+)(k,i))} . (33) 



is not effective due to the definition of the photon vacuum 
{Ai'^\k,t)) = (Sect. II). 



The equations of the second order for the correlators 
from Eq. ([33|) can be obtained easily with help of the 
equations of motion ((26l) - (f28l) . e.g., 



'd_ 

dt 



1 f Sk' ( 

+ i[uj{p^,t) + uj{p2,t) - k]Ub^{-p^,t)ap{p2,t)A\+\k,t)) =-ie(27r)-3/2y dV-^|^(p' - pi + k') 

\uvtp, (p', pi , k'; t) (a+ (p', t)ap (p2 , t)A.r> (k', i)A\+^ (k, t)) 
+ [vvtp' (P', Pi, k'; t){hp.{~p', t)ap{p2,t)Ar.{k', t)4+) (k, t)) 



S{P 



k!) 



[uuYp,p (P2, p', k'; t) (6„(-pi , t)ap, (p', t)Ar. (k', t)4+) (k, t)) 

H^'^(P2, P', k'; t)(fo„(-pi, t)6+ (-p', t)Ar, (k', t)AW (k, t))] } 
'5'a/3(pi,P2,k;t) + L/^^(pi,p2,k;t) . 



(35) 



On the r.h.s. of this equation there is a set of the terms which are absent in the standard QED without a strong 
originating from vacuum polarization effects in the pres- field. The group of terms 
ence of the quantized electromagnetic field (Sect. II), 
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C/^^(pi,P2,k;t) = C/^^^(p,,t)(6^,(_p,,t)a^(p2,t)4+)(k,t))-t/^y,(p2,0(&a(-Pi,i)a/?'(P2,^ 

+ UP^iput){a+{put)ap{p2,t)Ai+\k,t))~U^/Up2,t){bo.{~Pi,t)^^^ (36) 



can be omitted because of the approximation (|34p in or- 
der to close the chain of equations at this minimal level. 
The other set of terms leads to contributions of the vac- 

I 



The remaining processes can be easily identified in 
Eqs. ([35|) and ((37|) . where the random phase approxi- 
mation (RPA) is used as the truncation procedure. That 

I 



As the result, only the first and the fourth term in 
Eq. and the second term in Eq. ([57)1 survive in this 
order of perturbation theory (Approximation 2). 

The next approximation is the diagonalization of all 
one-particle correlation functions with respect to the mo- 
mentum variables and spin indices, 

{atip,t)ap{p',t)) = S^^Sip-p')f{p,t) , (40) 
(a„(p,i)4(p',i)> = (5„^<5(p-p'),r(p,i) , (41) 

These relations mean that spin effects are neglected. The 
analogous approximation was introduced for the photon 
correlation function, Eq. ([30]) . 

The processes of the instantaneous radiation of two 
photons have been omitted here, i.e., in Eq. ([35|l the 
substitution yl,,'(k',t) — > '(k^t) has been made (Ap- 
proximation 3). In order to rewrite Eq. (j35p in the inte- 
gral form, let us perform the intermediate transition to 



I 

uum polarization effects in the one-photon annihilation 
channel 



(37) 



leads to the following decoupling rule for the correlators 
of the type 



(38) 
(39) 

I 

the interaction representation 



a{p,t) — a(p,t)exp|i J dt'u}{p,t') j , 

b{p,t) = fo(p,i)exp|z^ dt'uj{p,t')Y (42) 
i(=^'(k,t) = A''^\k,t)exp{TiHt-ta)} , 



where io is some initial time. 

The approximations (f55)) - ([iT|) allow to rewrite the 
anomalous correlation functions from the l.h.s. of 
Eq. psp taking into account the vacuum polarization 



5;^(pi,P2,k;0 = -ze{27r)-'/'^l d'p\d^p'M-p'2-k) 

{[uu\a'i3' (P'l, P2. k; t){ba{-pi,t)af!{p2,t)a+ (pi, t)a^' (P2> *)> 
+ [uvjl,/^, (pi , P2 , k; t) (fe^ (-pi , <)a/3 (p2 , i )a+ (pi , t)bp (-Ps , i)) 
+ [^^*]a'/3' (pi : P2 ' k; t) (6„ (-pi , t)ap (p2 , t)6a' (-pi , t)ap, {p'2 , t)) 
+ Ma'/3'(pi> P2: k; t){bai-pi,t)ap{p2,t)ba'i-p'i,t)bp{-p'2,t))^ 

I 



{a+{p',t)ap{p2,t)Ar'{k',t)Al+\k,t)) ~ {a+{p',t)ap{p2,t)){Ar'{k',t)Al+\k,t)) , 
Ki-Pi,t)af,{p2,t)al'^\p[,t)b^^\~p'2,t)) ^ {be.i-p,t)b^pVip'2,t)){apip2,t)aiV{p[,t)) . 
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contribution (l37l) so that 



(6„(-Pi,0a^(P2,i)4+^(k,t)) = 
ie(5(p2 - pi + k) 



c^*'MQ/3(P2,Pi,k; t') 



x{[/(pi,t')+/(P2,t')-l][l+^r(k,0] 

+ [1 - /(pi,t')][l - /(P2,t')]}e-''^P^'P^'''^'''*' 
+C.C. , (43) 

where it was used that /'^ = 1 — / due to the electric 
charge neutrality of the vacuum ([1]) at t — > — oo and 

6'(pi,p2,k;i',t) = y" dr [a;(pi,T) + w(p2,t) - fc] . 

(44) 



The first group of terms in the curly brackets in Eq. ([43|) 
corresponds to the one-photon annihilation process (this 
contribution was investigated in the work [43|) while the 
second group describes the radiationless vacuum fluctu- 
ations. In the case of a strong subcritical laser field 
the number density of the radiated photons is small, 
Fr (k, t) <C 1 (Approximation 4) , so that the influence 
of the photon reservoir on the photon emissivity of the 
system can be neglected. Eq. (|43|) then takes the form 

(6,(-pi,t)a^(p2,t)AW(k,t)) - 
ie5{^2 - Pi + k) /■*^^,g-^e(pi,p„k;t'.t) 



2fc(27r)3/2 

X M^.^(P2, Pi, k; t')/(Pi, ^')/(P2, i') + c.c. . (45) 

Substituting ([^S]) into Eq. ([551) . obtain a closed ex- 
pression for the photon production rate 

X(p,p + k,k;t,t')/(P,t')/(P + k,t') 



-\-c.c. 



(46) 



where we have introduced the two-time convolution with 
respect to spin and polarization indices 



M;3a(P.P + k,k;t) 
[H;,3(P + k,p,k;t') . (47) 



is:(p,p + k, k;t,t') 



with the deflnitions (g]) and ((23)) . 

Additionally, it is assumed in Eq. (|46|) that the photons 
have equiprobable distributions regarding their polariza- 
tions, Fi= F2 = F. 

Thus, the photon production rate is a nonlinear 
(quadratic) non-Markovian function with respect to the 
electron-positron distribution function f{p,t). This or- 
der in the nonlinearity corresponds to the result [IJ] of 
the S-matrix approach. 

The consequent estimation procedure of the integrals 
on the l.h.s. of Eq. (|151) (method of photon count) was 
presented in [i^ on the basis of the methods given in [3! ■ 



The meaning of these approximations is the following. 
On the r.h.s. of Eq. (|46|) there is a high frequency mul- 
tiplier exp(— i6') with the phase ([M]). In order to select 
the low frequency component of the photon production 
rate ([46]) (this corresponds to the averaged, observable 
value), it is necessary to compensate this high frequency 
phase by means of the higher harmonics in the Fourier 
decompositions of the other functions in the integral (|46p . 
Indeed, the integrand in the expression (l44|l for the phase 
(mismatch) is very large 



w(p, r) + a;(p + k, t) — fc ^ 2m 



(48) 



and the energy conservation law is not fulfilled for the 
one-photon annihilation process. That is why the radia- 
tion of a real photon can be interpreted as a multipho- 
ton process. The photon number from the photon 
condensate of the external quasiclassical laser field with 
the frequency v can be estimated from the condition to 
compensate the mismatch (|48p by the energy of N^, quasi- 
classical photons. This leads to the estimate ^ 2m/iy. 
For optical lasers this is a huge number and therefore such 
kind of fluctuation event is very scarce. But in the case of 
7— radiation the result raises hopes for the possibility of 
an observation. In the general case, the eff'ect is defined 
by the external field parameters: the amplitude Eq and 
the frequency v. 

This conclusion about the role of multiphoton pro- 
cesses correlates with the analysis of the absorption co- 
efficient of the electron - positr on p lasma created from 
vacuum in the infrared region [4J|. According to the 
structure of the time dependent u, u-spinors [sj, the con- 
volution (|47p is a polynomial in eA{t) and hence it can 
not guarantee for the necessary compensation. Below it 
is assumed that the laser electric field A{t) = A^{t) = 
— [EQ/v)cos{vt) is subcritical, Eq < Ef.- Therefore, we 
use here the Markovian approximation (Approximation 
5) A'(p, p k, k; t, t') ^ K{p, p + k,k; t, t) = Kq = 4, 
see Appendix A. Thus there are two sources for compen- 
sation of the high-frequency phase (|44|) : the time depen- 
dence of the fermion distribution functions f{p,t) and 
the multiphoton contributions contained in the mismatch 

As it can be seen from the KE's ([T6|) . (fTTj) the fermion 
distribution function f{p,t) has two time scales defined 
by the frequency v of the external field and by the rest 
mass TO (or the one-particle quasi-energy uj{p,t)). Then 
the corresponding Fourier transformation can be defined 
by the double series 



f{p,t) 



fn.i{p)e 



(49) 



The numerical calculations have shown [20[ that the low- 
frequency behavior of the fermion distribution function 
in the presence of a linearly polarized field is defined ap- 
proximately by the second harmonic (Approximation 6) 



/(P,i) = /(P) [1 - cos(2z/t)] /2 



(50) 
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The higher low-frequency harmonics and the high- 
frequency harmonics in the decomposition ()49p are very 
small compared to ([SO]) and will be omitted in the fol- 
lowing. 

The choice of the following approximations depends on 
the value of the adiabaticity parameter [l^ 



7 = {EjE^)/{v/m) 



(51) 



The domain 7 3> 1 corresponds to the multiphoton pro- 
cesses in a rather small external field Eq <C Ec- In the 
case 7 <C 1 the external field is rather large, £'0 < Ec- 
These two limiting cases will be considered below. 



B. Multiphoton domain (7 3> 1) 

Since Eq Ec, one can use the appropriate perturba- 
tion theory in order to select the constant component in 
the photon production rate. The unique source is now 
the time dependence of the mismatch (|48p . For Eq <^ Ec, 
one can select in the quasi-energy the first order field ef- 
fect (Approximation 7) 



w(p,t) ~ wo(p) 



eEop 



cos(j/t) 



(52) 



i/cjo(p) 

where a;o(p) = uj{p,t)\A=o- Then the phase (|44|) will be 

e{t,t') c^no{p,'k){t-t') + a{p,k)[smiyt-smiyt'] , (53) 

where the mismatch in the absence of the external field 
is 



fJo(p, k) = wo(p) + wo(p + k) - fc 



and 



a(p,k) 



Eop , Eo(p-Kk) 



'^o(p) wo(p + k) 



(54) 



(55) 



The fast oscillating function in Eq. PS|) allows then the 
following representation 



cos$(t,i') ~ e 



iilait-t') 



^ J„(a)J„.(a)e'''("*-"'*')+c.c. , 



which is based on the known decomposition 

00 

exp(msiniy9) — J„(a)e' 



(56) 



(57) 



n— — 00 



where J„(a) is the Bessel function. Thus, the represen- 
tation (|56p contains two high frequency harmonics wo(p) 
and tJo(p + k) (sec Eq. ([54]) ). the set of the low frequency 
harmonics nv and the k- harmonic corresponding to the 
radiated photon. 

The following steps are these: we calculate the time 
integral in Eq. (|46p and select its time independent com- 
ponent corresponding to the observed value. Performing 



then the remaining momentum space integrations in the 
isotropic approximation and using the textbook formula 

,5[(/)(x)] =^{|0'(a;,)|}~'(5(x-a;,), ^{xi)^Q, (58) 

i 

we obtain then in the low frequency approximation 



= -rrr , ., ,„ ^«o+i 



(a) 



[Jno+sia) + Jno-i{a)] Kpi)f{pi + fc) , (59) 

where we took into account the lowest harmonics of the 
distribution function Eq. ((50|) and wi = ^/m^ + pf. The 
argument of the Bessel function in the isotropic approx- 
imation is 



2J¥^Eo 



Pi 



Pi 



Wl ^Luf + fc2 



(60) 



Let us explain the meaning of the momentum pi . By pq 
we denote the positive root of the equation f^o — = 0, 
where 



na{p,k) = ujo{p) + ^ujHp) + k'^ - k 



(61) 



is the mismatch and n = [fl(j{po, k)/u] {[x] is the integer 
part x) is the photon number necessary for overcoming 
the mismatch (l6T]). Then we obtain 



(nvf{nv-2kf , 
,{n.^kY 



1/2 



(62) 



Let now no be the minimal number of quasiclassical pho- 
tons, no ^ po = 0. The photon production rate ([62]) is 
equal to zero in this point. Let us suppose ni = no -I- 1 
in Eq. ([5^ . The momentum pi corresponds to this num- 
ber, pi = Til. Thus, according to Eq. (j63p we have 
/(fc) ^ Q!"°+^/^(0). It is difficult to estimate the general 
order of the perturbation theory, because the amplitude 
/(O) of the distribution function is calculated nonpertur- 
batively. 



1. The case of optical vacuum excitation (u <^m) 

In the optical part of the photon spectrum (fc < v) we 
have pi=\/km as the first root of the equation r2o(pi) — 
(no + 1)1' ~ (it corresponds to the leading contribu- 
tion from the set n > no), a = 2{Eo/ Ec){'m/iy)^/'^ < 1 
and no = [2m/v] ^ 1, i.e., the necessary number of qua- 
siclassical photons is huge and the intensity of photon 
radiation a very small. Indeed, from Eq. ()59p it follows 
that 



/(fc) 



1 dF(fc) 



m 
a 



dt 



Kr 





' Eq /m\3/2' 


2no "I 






J 



(.63) 
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The value /(O) can be estimated as a result of the numer- 
ical solution of the KE (|16p . (fT7|) for e~e~^ excitations in 
a laser field [45| . For the P W laser system Astra Gemini 
we have Eq - IQ-^Ec, X = 800 nm and /2(0) - 10"" 
[l^HOl- The spectral distribution of the radiated photons 
from the volume of the focal spot per time interval, 



dNk _ 
dtdk ~ TT^ 



I{k) k^m 



-m , 



(64) 



will be negligibly small for the mentioned parameters. 
However, the term in the curly brackets on the r.h.s. of 
Eq. ([63|) behaves as a 6'(a)-function with the branch point 
ao — 2 when letting a — ag. Then the spectral distri- 
bution ([64|) starts to grow strongly. Unfortunately, this 
value ag = 2 lies outside of the validity range of Eq. ([63| . 
Nevertheless, this gives a hint on the possible growth of 
the radiation intensity in this domain. Some additional 
analyses are necessary here. 

The 7-ray part of the photon spectrum k ^ m can not 
be considered in the framework of this approximation 
(a » 1 again) . 



2. The case of -y -ray vacuum excitation (u ^ m) 



In this case the mismatch (|61|) can be compensated by 
the smallness of the photon number from the quasiclas- 



sical laser field. 



no 



> 1 m 



Let no = 1 (this is the 



hypothetical limiting case^e.g., for the planned XFEL 
facility with A = 0.15 nm [IJI)- The developed theory is 
working well in this case (a ^ 1) for the subcritical fields 
E < Ec. 

In the optical part of the photon band (fc <C m) we 
have Pi = \/Zm and 

..v^|(^)\ (65, 

The spectral distribution according to Eq. (|59)) is 



dtdk 



2^ 



(66) 



Thus, the effect grows linearly with k. For v = \ MeV, 
Eo = IQ-^E^ and f{pi) - 10"" (see Fig. [J) we obtain 
again a negligibly small effect: the suppression factor is 
f{0)Eo/Ec ~ 10~^ so that a very weak signal results. 

The situation is slowly changing when going to higher 
frequencies of the excited signal (X-ray or 7-ray domain) 
at E/Ec = const. One can demonstrate this by writing 
pi ^ tor k^O and 710 = 1 

> , , .2 ^ 1/2 

! + k/m 



pi = m 



k/2m 



- 1 



(67) 



The situation becomes more optimistic at E ^ Ec 
when a ^ 1. For example, for the XFEL with E = 
0.24 Ec and A = 0.15 nm [ij] the intensity ([66]) can be 
accessible to observation, apparently. However, this case 
needs special investigation since for 7 1 the presented 
approach is not valid. 



C. Strong field case (7 < 1) 

We will use here the effective mass model [3l| based 
on the approximation 



uj{p,t) = y m'^ + {p ~ eA{t)) ->w,(p), 

with the effective mass defined by the relation 

ml = m"^ + ^ A^{t) + e'^E^/2iy'^ 

= 771^(1 + 1/27^), (69) 

where ^ ... >- denotes the time averaging operation, and 
7 is the adiabaticity parameter (|5ip . 

In this approximation the phase (I44[) becomes 
monochromatic 

0(pi,P2,fc;t',t) = r!4pi,p2,k)(t-i'), (70) 



174pi,p2,k) = tJ*(Pi) + w*(P2) - k, 



(71) 



i.e., the approximation (|68p leads to a suppression of mul- 
tiphoton processes (it corresponds to large values of the 
adiabacity parameter 7 3> 1) and the mismatch (I7ip can 
be compensated by the harmonics of the fermion distri- 
bution functions in Eq. (|46p only. 

The inspection of the fermion distribution function 
shows, in particular, that it oscillates basically with twice 
the laser frequency ([50|l . The substitution of Eqs. (f70|) 
and (|50p into the KE pS)) allows to perform the time 
integration, leading to the appearance of two harmonics 
in the radiation spectrum only (the 2"^ and the 4*''), 

F(k, t) = - (k) cos(2i/t) + (k) cos(4i/t) , (72) 



^(2)(k) = ^ / 7?^/(p)/(P + k)<5(2.-0.) 



2k 



(27r)3 



(73) 



(74) 

It is important that a constant component is absent 
here, because the mismatch ([IS]) could not be compen- 
sated in this case by other sources of the time dependence 
on the r.h.s. of Eq. (|46|) .^ 

Thus, in the case of the infinite system the solution ([72|) 
can be interpreted as "breathing" of the photon subsys- 
tem. However, the situation is changed, when the gener- 
ation of the e^e+7 plasma is considered in a small spatial 



^ This is in contrast to the case 7 S> 1, where accounting for multi- 
photon processes in the phase 1144 ll leads to a constant component 
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domain of the focal spot with volume ^ due to the vac- 
uum condition of the absence of the plasma in the 
initial moment of switching on the laser field. In this case 
one can expect, that all annihilation photons generated 
in the first half-period of the field will leave the volume 
of the system and therefore in the next half-period the 
reverse process (photon transformation to e~e~^ plasma) 
will be impossible. Such a mechanism leads to a pulsa- 
tion pattern for the photon radiation from the focal spot. 
It corresponds to introducing the condition of a positive 
definite photon production rate on the r.h.s. of Eq. ([72|l. 

For estimates of the amplitudes (|73p , ([7^ let us intro- 
duce the additional model approximation in the spirit of 
the model (l68ll. 



Lj^p + k) ^ uj,{p, k) = y/ujl{p) + fc2 , (75) 

and the isotropisation condition /(p + k) — >• f{p -\- k). 
The integrals on the r.h.s of Eqs. can then be 

calculated. For example, 

/l(2)^^ ^ f( \r( , ^ u}4po)uj4po, k) 
k cli*(po) + w4po,k) 

where 



As a representative characteristics of the effectiveness 
of the radiation from the focal spot domain we will con- 
sider the total photon number per time interval. 



dk k /(po)./(po + k)po . (80) 



The electron and positron distribution functions enter- 
ing here are defined as the solutions of the correspond- 
ing non-perturbative KE (Sect. II) describing vacuum 
creation of e~e^ pairs under the action of a strong, 
time dependent electric field of a standing wave of two 
counter-propagating laser beams. The cutoff parameter 
^max = 2to* is introduced in order to take into account 
the annihilation photons in the radiation spectrum. 

The fermion distribution function f{p,t) is a rapidly 
decreasing function with its maximum in the point p = 
[20I Isoj . On this basis for a rough estimate one can put 
Po = in the arguments of these functions on the r.h.s. 
of Eq. dlO]), 



^(,) ^ Mom. ^-^^^ dfcfc/(fc)po, (81) 



Po = 



i4iy^{k + i^y 



{k + 2vY 



(77) 



is the root of the equation — 
follows the threshold condition'^ 



2v = 0. From Eq. (|77l 



2v{k 



k + 2v 



(78) 



where according to Eq. ([77] 



Po(fc) 



9 

mi 



k + 4m* 



k 

56 

m* 



m / k 

— W48 + 56 

4 V m* 



15- 



fc2 



mf 



(82) 



This condition is rather nontrivial because the effective 
mass (|69p depends also on v. The minimal permissible 
value V = 2m* corresponds to fc = 0. For the 4**^ har- 
monic the threshold value falls to v = m^,, which is close 
to the parameters of the XFEL [13] . 

The 1/k - dependence on the r.h.s. of Eq. (|76|) cor- 
responds to the flicker-like noise of electrodynamical ori- 
gin. This feature in the spectrum of radiated annihilation 
photons has been found first in Ref. [i^ . 

The number of photons with the frequency k lying in 
the interval [fc, fc + dk] and radiated from the focal spot 
with the volume A'' = per time interval is defined by 
Eq. The fraction on the r.h.s. of Eq. ([7S|) is a slow 

function of the frequencies k and v and for the sake of 
a preliminary estimate it can be replaced by m*/2. For 
the 2"*^ harmonic we then obtain from Eqs. ([76]) and ([64]) 



dtdk 



/(Po)/(Po + k)po 



(79) 



gral ((8T|) . As the result, we obtain the following order of 
magnitude estimate 



(83) 



For the XFEL parameters Eq = 0.24 Ec and A = 15 
nm [1^ we have according to the kinetic theory in the 
e-e+ sector /(O) 10"^, see Fig.[Tl From Eq. JHS]) then 
follows 



10 



17 



(84) 



For the 4*^^ harmonic with the oscillation amplitude ([74]) 
the threshold for the generation of annihilation photons 
is lowered (see discussion after Eq. ([75]) ) but the intensity 
of the photon radiation is also lowered so that the order 
of magnitude of (|84p remains unchanged. 



SUMMARY 



^ A similar effect was found first in the tfieory describing the ab- 
sorption of a weak signal by the e~e+ plasma created from vac- 
uum [3. 



We have studied the photon production rate resulting 
from the one-photon annihilation mechanism in a quasi- 
particle EPF created from the vacuum under the infiu- 
ence of a strong electric field being a necessary condition 




FIG. 1: Amplitude of the photon distribution function /(O) 
as a function of the dimensionless laser field strength Eq/ E^- 

for the possible oecurrence of sueh a process [3l| . The re- 
quired strong electric fields occur, e.g., in the focal spot of 
two counter-propagating high-intensity laser beams. The 
methodic basis is an appropriately developed kinetic the- 
ory constructed in the quasiparticle representation. The 
fermion sector of the theory has been investigated be- 
fore [l^H^I on an essentially nonperturbative basis. The 
photon kinetics can be considered in the framework of 
the usual perturbation theory which allows to truncate 
the BBGKY hierarchy of KEs at the lowest order with 
respect to the fine-structure constant and to obtain the 
closed formula ([46|) for the photon production rate in- 
cluding vacuum polarization effects. 

We have investigated this expression for different char- 
acteristics of the laser field (frequencies v and amplitudes 
Eq) which is conveniently described with the adiabaticity 
parameter 7. In order to become observable it is neces- 
sary for the photon radiation process to compensate the 
energy mismatch (|48|l by a sufficient number of quasi- 
classical photons from photon reservoir of the external 
electric field. The strength of the laser field Eq defines 
the intensity of this process. Thus, it can be understood 
that in the case of optical lasers the number of quasiclas- 



sical photons in the optical range needed for the compen- 
sation of the mismatch (|48|) is huge and the intensity of 
such process is negligibly small. 

However, the intensity of photon production increases 
strongly in the domain of 7— ray excitations of the vac- 
uum (where only a small number of quasiclassical pho- 
tons is necessary to overcome the mismatch) and reaches 
the values sufficient for experimental observation, e.g., 
for the projected XFEL [ij]. This result is qualitatively 
confirmed by the recent findings of Ref. [1^ which have 
been obtained within a different approach. 

In the case of the two-photon annihilation channel the 
necessity to compensate the mismatch of the type (|48|) 
is absent and therefore the photon production rate can 
be large enough to be observable already in the optical 
domain. This channel of the EPP annihilation will be 
investigated separately. 

We would like to note also that the type of kinetic the- 
ory based on the quasiparticle representation, which has 
been described in the present work can also be useful for 
the investigation of other reaction channels in the e~e"'"7 
plasma. One of such channels is the cascade process of 
EPP multiplication in a strong electric field (e.g., [soj). 
The example of one-photon annihilation presented here 
shows, that strong external fields can lead to qualitative 
modifications of processes relative to standard QED. 
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Appendix A: Convolution (|47p 

As an example, we will consider here the calculation 
of the convolution ()47p in the Markovian approximation. 
Setting t = t' , we obtain 

K{pi,p2,k;t,t) = 

(Al) 



With help of the relations [5l| 

Y^[Uaip)]k[Ua{p)]i = 

a 



1 



Mp) 
1 



Mp) 

one can transform Eq. (jAip to the form 



[p + m]kt, 
[p-m]ki 



(A2) 



K{pi,p2,'k;t,t) = 

^ -Tr(pi - m)r{p2 + m)YeUk)el{k) 



4lLj(pi)w(p2 
1 

UJ(PI)UJ{P2) 



{2(pie'-)(p2e^) - {m^+pw2) {e^elg'^'')] 



(A3) 



The photon polarization vector e^(k) is satisfying to the 
conditions {r — 1,2) 



nv r r 

9 e^e^ 



-2, 



(Pie'')(p2e'') = (pie'')(p2e'") . 



(A4) 
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Thus, we finally obtain 



Since in the integral we have 



A'(pi,p2,k;t,t) = 
2 

{w? +P1P2 + (pie'')(p2e'')} .(A5) 



In the case pi = P2 = it follows from (|A5 
ii:(pi,P2,k;t,t) = A'o ==4 . 



(A6) 



The last formula is used in the qualitative estimates of 
the photon production rate in the different models. 
For Pi ^ P2 7^ one can obtain another estimate. 



P1P2 => 3'5»fcPiP2, 



we obtain from Eq. (jA6l 



(Pi,P2,k;i,t) 



w(pi)w(p2 



+ uj{pi)uj{p2) 



-P1P2 } .(A7) 



